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COMMENTS ON DIFFERENCE SCHEMES FOR THE THRFE-DIMENSIONAL 
TRANSONIC SMALL-DISTURBANCE EQUATION FOR SWEPT WINGS 

By .Jerry C. South, Jr. 
l.anqley Research Center 

SUMMARY 

Certain problems arise In constructing stable finite-difference schemes 
for the three-dimensional transonic smal 1 -disturbance equation when cross- 
flow terms are Included to better »pproximate swept wings. These problems 
are discussed and some possible remedies are offered. 

INTRODUCTION 

Recently, some effort (refs. 1 to 3) has been devoted to improving the 
three-dimensional transonic small disturbance equation for applications to 
swept wings. It seems that the Important crossflow term that needs to be 
Included Is <|) , where x 1s the longitudinal coordinate and y Is the span- 

y 

wise coordinate. Without this term, the small disturbance equation cannot 
adequately predict weak, swept shock waves. The purpose of this paper Is to 
point out some numerical difficulties that may arise In constructing stable 
finite-difference approximations to equations which include this additional 
term. 
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ANALYSIS 

( Jameson (ref. 4) has sho<n that it is advantaqeous to examine the 
canonical form for the potential equation, particularly when studying 
appropriate methods for combining central and backward differences for various 
derivatives in supersonic regions of the flow. As in reference 4, then, the 
canonical form is written as: 


where 

ss 


(a^-q*) + a* - 0 

represents differencing in the local streamline direction, and 

q" - 1 + 24.^ + -|)J t 4,2 + 4,2 

^ (q*'i) 

a* * M'* 

00 00 

q*4i ■ u*4> + v*4» + w*4> 

^ ^ss ^xx ^yy ^zz 

+ ?uv4)^y + 2 uw4>j^, + 2vw4)y^ 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


u = 1 + 4> 

X 

(6) 

V • <Py 

(7) 

w = 4>^ 

(8) 
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Combining equations (2) to (4) gives: 


a* - q* 


1-M» - (q* - 1) 

m 2 nr» * 


' - -LJii m: (q’ - ' ) 




If terms Involving perturbation velocity products are Ignored, equations 
(2), (5), (9), and (10) become: 


q* ^ 1 + 24 ., 


a* - q' 


1 - - (y+1)M^4) 

00 ' oo^x 

CO 


V . 1 - (y-1 )m’* 


(b ^4 +2d>4 +24>4> 

^ss *^xx y xy *^z*^xr 


The term V*4) - is : 


7^4 “ d> * (b ■^4 " 24 4 “ 24 4 

^ ^ss ^yy ^zz y^xy *^z^x2 


Substituting equations (12) to (15) into (1) yields: 


A (4 + 24 4 + 24 4 ) 

XX y xy ^z*^xz 


+ B (4 +4 - 24 4 - 24 4 1 » 0 

'^yy ^zz *^y\y ^z^xz' ^ 
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•Vre 


(17) 

(18) 

In the Interest of further simpl icatlon, nnst researchers have Ignored 
the term ** unimportant for the swept wing problem. Some discussion 

of the relative size of this term Is given In reference 2. For the purpose of 
this paper, deleted from equation (16), yielding finally: 

^ ^ ^ ^zz ■ ^Vxy^ * ° 

where second derivative terms belonging to are barred to Indicate that 
they shoulo be represented by upwind differences In supersonic regions of 
flow, a'.id terms deriving from 7*ij) - should be represented by central 
differencer. In the spirit of Jameson's rotated difference scheme (ref. 4). 

To a close approximation, equation (19) changes from elllpt'. ■ hyperbolic 
type when the coefficient A changes from positive to negative. Hence, a 
stable numerical scheme with artificial viscosity of the correct sign Is 
obtained by using upwind differences for ^7“ and whenever A<0- 
That Is, for we have: 


s - 1-M» - (y + 1) M'<^ 
8 - 1- (y - 1) 


Ax' 


XX 


•^l.j.k ■ 2^1-1, j,k ^ <^1-2,j,k 


( 20 ) 
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and for we have: 

A/ 


If ♦y > 0 


and 


AxAy s^y ^ - ♦i.i.j.ic + ♦i-l,j-l,k 


( 21 ) 


AxAy (>j^y - <>i j+i,k • '^l-lj+l.k ♦l-l.j.k 


- <j) 


U.k 


( 22 ) 


if ♦y < 0 

It 1s Important to note that such formulas always enhance diagonal dominance; 
that Is, they produce a contribution to the diagonal (coefficient of ♦iii.) of: 


A 



+ 



(23) 


Such schemes are highly desirable from the viewpoint of numerical stability. 

Equation (19) still has an important deficiency: because of the absence 

of certain deleted small terms, it cannot be cast into divergence form. 

Further expedient approximations can be made, however, to achieve a divergence 

form. For example, in reference 1, the total coefficient of was 

y ^y 

collected and approximated by ignoring 


2(A - B) « (I + 2^^) ^ -ZM^ 


(24) 
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and the - contribution In the coefficient of ♦ was also Ignored, yielding: 

A i - ?M* AA+Bd>+tfi"0 (25) 

^xx • ^y^xy ^yy ^zz ' ' 

The bar over ^ has been dropped In equation (25) because of the slight 
xy 

confusion created by combining the coefficients of and In equation (24). 

xy xy 

strictly speaking, to retain the Identity of and ♦ In equation (25), one 

xy xy 

would write: 


A 4 + 2(1-M*) t <ti -2(t> (> + + (1> ■ 0 (26) 

^xx ' <"' y xy ^y xy yy zz ' * 

However, upwind differencing for , as given by equations (21) and (22), 

xy 

detracts from diagonal domlnancel That Is, application of formulas (20) to 
(22) In equation (26) produces a diagonal contribution of: 


0-mM I a 


Ax’ 


Ax Ay 


(27) 


and the two terms compete with one another, since A<0. For large cross flows 

(large |l>y|), such a procedure would be unstable. From a numerical point 

of view, then, one should use an upwind difference only for and use 

central differences for all other terms In equations (25) or (26). An 

alternative is to retard In the x-direction only, for example: 

xy 


2AxAy«^y = *^i,j+l,k ■ ♦i-l.j+l.k ^ '*’i-1.j-1.k 


but then that is not a ‘‘rotated'* difference schenip, since the direction of 
the crossflow is ignored. Furthernv're, the artificial viscosity that is intro- 
diced has the wrong sign, since its coefficient is 1 - rather than A. 

In reference 2, this confusion was avoided by using a different 
approximation to equa^^on (19). Namely, the term ^ is dropped, and B 
approximated by 1 , so that equation (19) becomes: 

" ♦yy ‘ ♦zz • Vxy " “ («) 

Here there is no question about the differencing of the term, since it 

xy 

may be considered as deriving entirely from thus would always 

be centrally differenced. 

In reference 3, the difficulty of constructing a divergence form after 
neglecting certain terms is overcome by a simple but effective idea: the 
equation is expanded in the original conservative form, 

[p(l+ii>^)]^ ■ 0 (30) 

where p is the density. The terms neglected in reference 3 are essentially 

the same as those dropped in reference 1, and the same delicate point finally 

arises concerning and 0 . It is stated that retardation in the x-direc- 

xy xy 

tion only is then used for the ({(^^-contributions, so a formula related to 
equation (28) is used. The equivalent artificial viscosity will have the wrong 
sign, for the same reason as before. Although the title of reference 3 Includes 
the word "rotated difference," the claim here is that this is not a rotated 
difference scheme, but rather a "split" scheme, as in reference 2. 
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CONCLUDING RFMARKS 


Experience has shown that rotated difference schemes have nice stability 
properties, since If they are carr’^d out correctly, each upwind difference 
(Including cross derivatives like 4 . ) enhances diagonal dominance, as 

Illustrated In equations (21) to (23). This feature Is an Indication that the 
artificial viscosity Is correct Iref. 4). Retarding a difference formul* for 
♦ in the x-directlon only does not Insure the correct sign to the resulting 

Aj 

artificial viscosity unless the coefficient of 4> Is the correct sign. Such 

/Cjr 

retardation might damage stability more than a central difference formula, 
which produces zero viscosity. 

It Is felt that the Dutch apprt'ach (ref. 3) Is a good one for deriving 

a conservative approximation for swept wings, but that more terms should be 

retained to force a "switching" coefficient similar to A (eg. (17)) on the 

equivalent ^ . In this way, a rotated scheme can be devised with correct 
xy 

art^flcal viscosity, and the concomitant diagonal dominance. Otherwise one 
should use a central difference formula for 

xy 
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